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ABSTRACT 

The recent discovery that the Fe-K line luminosities and energy centroids observed in nearby SNRs 
are a strong discriminant of both progenitor type and circumstellar environment has implications for 
our understanding of supernova progenitor evolution. Using models for the chemical composition of 
core-collapse supernova ejecta, we model the dynamics and thermal X-ray emission from shocked 
ejecta and circumstellar material, modeled as an r~ 2 wind, to ages of 3000 years. We compare the 
X-ray spectra expected from these models to observations made with the Suzaku satellite. We also 
model the dynamics and X-ray emission from Type la progenitor models. We find a clear distinction 
in Fe-K line energy centroid between core-collapse and Type la models. The core-collapse supernova 
models predict higher Fe-K line centroid energies than the Type la models, in agreement with obser¬ 
vations. We argue that the higher line centroids are a consequence of the increased densities found 
in the circumstellar environment created by the expansion of the slow-moving wind from the massive 
progenitors. 

Subject headings: ISM: supernova remnants — stars: mass-loss — supernovae: general — X-rays: ISM 
— ISM: abundances 


1. INTRODUCTION 

Supernovae, both thermonuclear (la) and core¬ 
collapse, represent the endpoints in the stellar evolution 
of either Chandrasekhar white dwarfs or massive stars 
with zero age main sequence masses > 8M 0 . Core¬ 
collapse supernovae (CCSNe) exhibit a wide diversity 
of properties. Classifications are based on spectra and 


nism (Filippenko 

1997) 

by the presence (1 

1) or 


trum. Type I CCSNe include Type Ib/c, Ibn, and Ic- 
BL. The progenitors of these supernovae are thought 
to be either mas sive He stars or som e subset of Wolf- 
Rayet (WR) (see Gal-Yam et al.|[2007 for a broad clas¬ 


sification). Type II SINe include the subtypes IIP (e.g.. 
SN 1999em: ISmartt et al.||2002l) IIL (e.g., SN 1979C: 


Fesen fc Matonick|| 1993 1, ll-pec (e.g., SN 1987 A; |Cheva- 


fier fe Emmcring 1989), Iln (e.g., SN 1995N; Fransson 


et al. 2002), Hb (e.g., SN 1993J; Fransson et a T 1990), 


and possibly some superluminous supernovae that ex- 
hibit evidence for a strong c ircumstellar interaction, such 
as SN 2006tf ( [Smith et al.||2008 >. 
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Unlike Type I CCSNe, the progenitors of Type II su¬ 
pernovae are likely red or yellow supergiants (R/YSG), 
or in the case of Iln, luminous blue variables (LBVs; 
Gal-Yam et al. 20071. In at least a few cases, progeni¬ 
tors h ave been identified in pre-explosion images ( Smartt] 
2009). In both Type I and II CCSNe, the progenitor 
mass-loss rate can vary between 10~ 6 10 -4 M 0 yr _1 

(and for Ilns and SLSNe, 0.1- 10 M 0 yr _1 , presumably 
as eruptive mass-loss), but the wind velocity ( v w ) be¬ 
tween the RSG or WR progenitor can differ by as much 
as two orders of magnitude ( v w = 10-20 km s - 1 in a 
RSG, and ~ 1000 km s” 1 in a WR star; Crowther]|2007| 
Smith 2014). For steady, isotropic mass-loss, the circum- 


stellar medium (CSM) density is proportional to M/v w . 
Thus high mass-loss rates coupled with slow wind veloc¬ 
ities can lead to a substantial amount of mass close to 
the supernova progenitor. This would be in contrast to 
the large, evacuated cavities expected around WR pro- 
genit ors, du e to thei r large mechanical wind luminosities 
(c.f., |Koo fc McKeefl992] ) . 

While typing supern ovae based on th eir optical spec¬ 
tra is well established (Filippenko 1997), connecting su¬ 
pernovae types to supernova remnants (SNRs) remains 
difficult. In the case of Type la supernovae, their rem¬ 
nants can frequently be readily identified by the iron 
content observed in the X-ray spectrum (|H ughes et al. 


1995 Badenes et al. 2006| 2008 Pat naude et al.||2012[ ). 


Given the broad diversity of CCSiN types, connecting 
remnants to supernovae (and to supernova models) is 
a challenge. However, examples exist where supernovae 
have been linked to supernova remnants, or supernova 
remnants have been directly typed. For instance, light 
echo analysis allows for direct typing of some supernova 
remnants such as Cas A ( Fransson et al.|1996[ Krause et 


al.|2008||Rest et al.|2008|), andlMilisavljevic et al.|(|2U12|) 

recently showed that the [O IIIJ line shape in SN 1993J 
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was remarkably similar to that of Cas A. Finally, simi¬ 
larities in X-ray spectra provide evidence that SN 1996cr 
ma y be a Type II pec, similar to SN 1987A (Bauer et 
aL|2008l. However, beyond these examples, connections 
between SNe and SNRs remai n sparse, particular ly for 
the class of core-collapse SNe (Hughes et al.|1995|. 


Recently, Yamaguchi et al. (2014) presented a method 


of typing supernova remnants based on the Fe-Ka line 
centroid and luminosity. Since Fe is produced in the 
center of the progenitor during the explosion, heating of 
Fe can be delayed, resulting in an ionization state lower 
than He-like (Fe 24+ ) in young and middle-aged SNRs. 
The ionization state affects the Fe-K line centroid, which 
can be m easured to high precisio n with current X-ray 


satellites. Yamaguchi et al. 


empl oy the de la yed d etonation models DDTa and DDTg 
(Badenes et al. 2003). In these models, a flame propa¬ 
gates as a slow deflagration, with a transition to a det¬ 
onation induced at some flame density. The transition 
point sets the amount of 56 Ni produced in the explosion. 
The DDTa model produces ~ 1 M 0 of 56 Ni, while the 
less energetic DDTg model produces only ~ 0.3 M 0 of 
56 Ni. 

To model the core-collapse supernovae, we employ a 
range of ejecta models from a variety of single star pro¬ 
genitor scenarios. For this paper we have computed new 
stellar models of 12 M 0 and 25 M 0 initial mass, Mod¬ 
els sl2D and s25D. These models have b een computed 
using the KEPLER stellar evolution code (Weaver et al. 


(20141 showed that the Fe-K 19781 and are similar to those in Woosley fc Heger 


use th e updated solar abundances from |Lodders et al.| 
(2009). For both models, the stellar evolution, inclucU 


line centroids tor la SNRs are generally lower (< 6550 (see also Woosley et al.|[2002 Rauscher et al. 
eV) than those found in core-collapse SNRs. Addition¬ 
ally, they found that when computing synthetic Fe-K line 
centroids and luminosities from well tested mod els f or 
Type la ejecta, the models fBadenes et al.||2006[ |2 008|) 
predicted bulk properties in line with observations (j Ya- 
maguchi et al. 20141. 


2002), but 


Here we extend the work of 
address the question of whet' 


Yamaguchi et al. (2014) to 
rer models for UCSiNe are 


able to reproduce, in broad terms, the observable bulk 
properties of core-collapse SNRs. We employ a model 
which tracks the hydrodynamics and time-dependent 
ionization of shocked circumstellar material and super¬ 
nova ejecta, coupled to an emissivity code to compute 
Fe-K line centroids and luminosity, as a function of super¬ 
nova ejecta model, age, and circumstellar environment. 
We compare our models to Suzaku observations of Galac¬ 
tic and Magellanic Cloud SNRs, and discuss implications 
and future directions which this work can take. 

2. HYDRODYNAMICAL MODEL 

We employ our CR-hydro-NEI code which allows us 
to simultaneously describe the thermal and non-thermal 
emission at the forward and reverse shocks in young 
SNRs. CR-hydro-NEI is a ID Lagrangian hydrodynam¬ 
ics code based on VH-1, a multidimensional hydrody- 
n amics c o de developed b y J. Blondin and colleagues (e.g., 
Blondin & Lufkin 1993). CR-hydro-NEI simultaneously 
models the supernova blastwave dynamics and particle 
acceleration, including the back reaction of the nonther- 
mal particles on the SNR dynamics. It is capable of 
modeling the nonthermal particle spectrum as well as 
the broadband nonthermal and thermal emission from 
shocked circumstellar material and ejecta. Here, we run 
the model without considering the effects of diffusive 
shock acceleration on the dynamics and emitted spec¬ 
tra (i.e., we set the particle injection to the test particle 
limit). S pecific details concern i ng CR-hydro-NEI can b e 


ing mass-loss, is followed to the point of core-collapse. 
Model sl2D loses approximately 3 M 0 of material, while 
Model s25D loses ~ 8 M 0 due to stellar winds by the 
time of core-collapse. Details of the explosion can be 
found in Rauscher et al. (2002f). As a brief summary, a 
piston is placed a the base of the oxygen burning layer 
where entropy per baryon rises above S/ks = 4. The 
piston is first moved inward at about 1/4 of local gravi¬ 
tational acceleration until a radius of 100 km is reached; 
then it is moved outward at a constant fraction of local 
acceleration until a radius of 10,000 km is reached. The 
outward acceleration is adjusted such that an explosion 
energy of 1.2xl0 51 erg is reached; the estimated explo¬ 
sion energy for SN 1987A is taken as a typical value here. 
T he prescription for mixin g in the ejecta can be found 


Heger & Woosley (2010), but the default parameters 


used here have been adjusted such that enough 56 Ni is 
mixed into the envelope to reproduce the light curves of 
Type II supernovae, in good a greemen t with hydrody¬ 
namic simulations (Joggerst et al. 2009). 

We also use a model tor ( JCSN ejecta tailored to 


SN 1987A (Saio et al. 1988 

Hashimoto et al. 

1989 

Shigeyama &; Nomoto 

1990). 3 

L'he model comprises a 6 


Finally, we include a model that has been previously ap - 
plied to the Type lib SN 1993J ( Shigeyama et al.|1994 ). 


found in Elli s on et al. ( 20071; iP atnaude et al. ( 2009); 
Elli son et al.l (|201U|); IFatnaude et al.|(|2010D; and|Lee et 
al.l (I2014D. Recent modifications to the code al 


The progenitor for this model is an 18 M 0 mam sequence 
progenitor with a metallicity Z = O.O2Z 0 . This progen¬ 
itor loses « 15 M 0 of material prior to core-collapse. 
Further d etails of the mode l s are s ummarized in Table [lj 
Though Yamaguchi et al. (2014|) have previously shown 
agreement between Type la models and measured Fe-K 
line centroids and luminosities, we include a subset of la 
models in our sample to validate our results. For rea¬ 
sons discussed below, we do not expect a n exact agr ee- 
ment between our la results and those in Yamaguchi et 


ow us 

to track the evolution of the thermal emission from the 
shocke d ejecta, as well as employ custom ejecta models 
(Lee et al.| |2014 |. 


al. (2014), but broad agreement is provided as a con- 


2.1. Ejecta Models 

For this study, we make use of previously available 
models for the composition of the ejecta in both core¬ 
collapse (CC) and Type la SNe. For the la models, we 


sistency check. The core-collapse SN ejecta models are 
from the evolution of single-star progenitors and thus do 
not include the effects that binarity have on the progen- 
itor e volution and the circumstellar environment (Sana 
et al.| |2012). We will address the role of binarity on both 
the progenitor and CSM in a future paper. 

2.2. Circumstellar Environments 
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A key question is the shape of the circumstellar en¬ 
vironment around the supernova progenitor. Red and 
yellow supergiants (R/YSG) expel several solar masses 
of material over their lifetimes (10 -6 10 -4 M 0 yr _1 ) 

with velocities of 10-100 km s _1 . Wolf-Rayet progeni¬ 
tors expel similar amounts of mass per year, but with 
wind velo cities 100x that of their R/YSG counterparts 
(e.g., S mith|20~T4 ). Additionally, episodic mass-loss may 
be re levant in the th ermal and dynamical evolution of 
SNe (Pan et al.|2013 ), resulting in emission that is lumi¬ 
nous early on, but decays sharply due to a sudden de- 


CSM and ejecta as a function of SNR age. We pass these 
parameters to a n emissivity code (Patnaude et al.|[2010 


crease in the circumstellar de nsity (e.g., Dwarkadas 2011, 
Dwarkadas & Gruszko 2012). Regardless of details, it is 


currently accepted that SNe of the Type II variety expel 
material with velocities ~ 10 - a few hundred km s _1 , 
while core-collapse of the I b/c (a nd variants) have wind 
velocities ~ 1000 km s _1 (Smith 2014). 

Given the diverse nature of the CSM environment ex¬ 
pected around CCSNe progenitors, we chose to take a 
general approach and consider a range of mass-loss rates 
and wind velocities appropriate for massive, red super¬ 
giant progenitors. We consider progenitor winds with 
mass- loss rates of M = 1 - 2xl0 -5 M 0 yr _1 , and wind 
velocities of v w = 10-20 km s _1 . We do not consider 
CSM environments shaped by fast winds, e.g., Wolf- 
Rayet and O-type stars, as these win ds are expected t o 
result in low density cavities (see, e.g., Dwarkad as|2 007). 
The low density CSM does not lend itself to substantial 
thermal X-ray emi ssion (e.g., as in SNR RX J1713.7-3946 
Ellison et al.|2012 ), and any shocked CSM has a low tem- 
perature due to the long Coulomb heating timescale. For 
the la models, we assume uniform densities with ranges 
of 0.3 - 3.0 cm -3 , consistent with results inferred from 
observations o f t he e nvironments around la progenitors 
( Badenes et al.|[2007 ). 


3. MODELING AND RESULTS 
We evolve the C C SNe to ages of 3000 years 


Dwarkadas & Chevalier (1998) pointed out that for a 
Type la ejecta model interacting with a constant density 
ambient medium, the reverse shock propagates to the 
center of the ejecta when the forward shock has swept 
up ~ 24x the ejecta mass. For the la models consid¬ 
ered here, the reverse shock reaches the ejecta center at 
ages of ~ 1000 years for the densest CSM environments. 
Thus we only consider the evolution of the Type la mo d¬ 
els to 1000 years. This differs from |Yamaguchi et ah] 
(2014) where they evolve the la models to much larger 


ages, allowing the shock to bounce and re-shock previ¬ 
ously shocked ejecta. For both the CC and la models, we 
evolve them over the range of parameters summarized in 
Table [l] resulting in a grid of models that span a range 
of circumstellar environments. 

The ejecta models include not only the composition 
as a function of mass coordinate, but also the density 
and velocity as a function of mass coordinate. For this 
study, we have mapped the ejecta composition onto ap¬ 
proximations of the ejecta structure: the core-collapse 
supernova remnant models assume ejecta with a co nstan t 
den si ty core and an n = 9 power law envelope (Truelove 
& McKee|1999). The Typ e la ejecta are modeled with an 
exponential ejecta profile ("Dwarkadas & Chevalier 1998). 
Our code produces ionization tractions tor both shocked 


Lee et ak|2014) to compute the fu ll SNR spe ctrum as a 


function of time. As discussed in Lee et al. (2014), we 
include thermal and Doppler broadening, though we do 
not consider their effects on our results here, nor do we 
consider the individual contributions from shocked CSM 
versus shocked ejecta as a function of time. Future X-ray 
observatories such as Astro-H may be able to discern be¬ 
tween these two components. At chosen time-steps, we 
used Xspe<j/]to synthesize a Chandra ACIS-S observation 
of the SNR. We then fit the 6 - 7 keV emission from the 
simulated observation to a Gaussian with a power-law 
continuumj^] The results from these simulations, for the 
grid of Type la and core-collapse SNe models lists in Ta¬ 
ble [l] are shown in Figure l] where we plot the Fe-K line 
centroid as a function of luminosity (left panel), as well 
as the evolution of the Fe-K line luminosity as a function 
of scaled SNR age (right panel). 

4. DISCUSSION 


The data i n Figure |T| (left) are taken from Yamaguchi| 
et al. (2014|. In broad terms, we find results that are 


consistent to those in Yamaguchi et al. (2014): super¬ 
nova remnants typed by their Fe content as las are gen¬ 
erally consistent with the synthesized X-ray spectra from 
la models. Note that we do no t expect an exact match 
between the results presented in Yamaguchi et al. (20141 
and here: our calculations include emission from both 
shocked ejecta and circumstellar material, while they 
only consider emission from shocked ejecta; the electron 
to proton temperature at the shock in their model is as¬ 
sumed to be either 0.01 or 0.03, while we assume mass 
proportional heating at the shock, follo wed by Co ulomb 
heati ng downstream; as discussed in (Patnaude et al. 
2010), the collisional ionization and recombination rates 
used in the ionization balance calculations, as well as the 
atomic data used in the spectral synthesis likely also dif¬ 
fer. However, even with these numerous differences, it is 
encouraging to fin d agreement between our la results and 
those presented in Yamaguchi et al. (2014), and gives us 
confidence that the results for the core-collapse models 
are insensitive to these differences. The line centroids 
and luminosities for the core-collapse models are also 
presented in Figure Jl] (left). For those models, we find 
that the measured Fe-K luminosities and line centroids 
are consistent with the models, with a few notable ex¬ 
ceptions discussed b elow . 


As discussed in Yamaguchi et al. (2014), all la SNRs 
have observed Fe-K line centroids < 6550 eV, and the 
chosen la ejecta models also do not produce Fe-K cen¬ 
troids in excess of 6550 eV. Even Kepler’s SNR, which is 
thoug ht to be the result of a luminous la in a mo dified 
CSM (Chiotellis et al. 2012 Patnaude et al. 20121 still 
displays bulk properties consistent with a la origin. 

For the observed core-collapse SNRs, the observed 
line centroids are broadly consistent with the model pre¬ 
dicted centroids and luminosities. In most cases the core¬ 
collapse models show Fe line centroids that are in excess 

7 http://heasarc.gsfc.nasa.gov/xanadu/xspec/ 

8 Since we are fitting the simulated spectrum over such a narrow 
band pass, the continuum, which in reality is thermal in nature 
and the sum of several shocked components, is well approximated 
by a power-law. 
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of 6550 eV. Only in the youngest models (fsNR < 200 
yr) d o we fi nd Fe line centroids at or below 6550 eV. 
Yamaguchi et al. (2014) have argued that the higher ob¬ 
served line centroids are likely the result of higher am¬ 
bient medium density, and our results support this con¬ 
clusion. For the highest density la model, n arn i, = 3.0 
cm -3 , while for the core-collapse models, even though 
the density follows a r ~ 2 power law, n am b > 10 cm -3 at 
radii in excess of 10 17 cm, owing to the high density close 
in to the progenitor. The higher ambient medium den¬ 
sity has a two-fold effect: first, the shocked CSM is at a 
higher temperature and ionization state, since heating oc 
n 2 , and ionization rates are oc n; secondly, the dense cir- 
cumstellar environment produces a strong reverse shock 
which can deposit more energy into the ejecta, resulting 
in higher overall charge states in the shocked ejecta. 

In Figure [I] (right), we plot the time evolution of the 
Fe-K emission from both the forward and reverse shocks, 
for the range of mass-loss parameters listed in Table [l] 
W e scale the SNR age to the c haracteristic SNR age given 
by Truelove & McKee (19991: 


to the progenitor. Outside of this is a low density bubble 
created by the fast wind from the main-sequence phase. 
Models s!2D and s25D r epresent such RSG progenitors 
(Heger & Woosley||2010|). Model s25D has lost almost 
: 3iVl 0 of material over the course of its evolution, indi¬ 
cating that the RSG phase for this progenitor cannot be 
more than ~ 10 6 yr, though it is probably less than this. 
The extent of the region occupied by the RSG wind is set 
by the pressure of the surrounding interstellar medium, 
and can range from < 1 pc in IIP proge nitors to greater 


than 5 pc in Ilb/L progenitors (|Chevalier 


2005). The 


U = 17702^ 


3/2 


M. 


w \v w ,6 yr , 


(1) 


where the ejecta masses and explosion energetics are 
taken from Table [lj As seen in Fig. [l] (right), the emis¬ 
sion from swept up circumstellar ana shocked ejecta is 
initially comparable in the sl2D, s25D, and SN 1993J 
models, but the emission from shocked ejecta quickly 
surpasses that of the shocked circumstellar material by 
a factor of two or more. Interestingly, and previously 
pointed out in Lee et al. (2014), there is no apprecia- 


modeled radii shown in Figure [2] are generally consistent 
with this, and may suggest that at radii in excess of 10 
pc, the wind from an earlier phase of evolution is required 
- as discussed in Section 2, we applied these models with 
a broad parameter space, and more detailed modeling for 
each individual object is necessary. 

4.1. Outliers 

As seen in Figure [I] (left), while there is generally good 
agreement between the supernova explosion models and 
the bulk properties, several core-collapse measurements 
stand out as having higher than average Fe-K line lumi¬ 
nosities: Cas A, N132D, W49B, and N63A. In the case 
of Cas A, it appears likely that there was bulk overturn- 
ing of the deepest layers of ejecta during the explosion 
(Hughes et al. 2000). This resulted in the Fe-rich ejecta 
being brought to the surface during the explosion. Sub¬ 
sequently, the iron is observed to have been shocked at 
early ti mes, resulting in a high charg e state and lumi¬ 


nosity (Hwang & Laming 2003, 
SNRs M l 321) and JN63A are both 


2012). The two LMC 
thought to be interact¬ 


ing with a large amount of interstellar material in star 


models, even at late times (our models do show that the 

material l 

Warren et al. 

2003 

Borkowski et al. 

2007 

)• 

Fe-K emission from shocked ejecta does begin to appear 

W49B also displays 

ngher than average Fe-K lumi- 


at late times, though remains very weak). The SN 1987A 
model differs from the other models in two aspects: first, 
it is a compact blue progenitor, and the densities in the 
ejecta are higher than in the other models with extended 
H-rich envelopes, and secondly, the Fe abundance in the 
SN 1987A model is ~ two orders of magnitude lower than 
in the other models- there is not a significant amount of 
Fe to shock in the ejecta. 

We also investigate the relationship between the Fe-K 
line luminosity and energy centroid and the SNR radius. 
In Figure [2] we plot the line luminosity as a function 
of radius for each model (left panels) while in the right 
panels we plot the line energy centroid as a function of 
radius. Note that for those models with a high mass- 
loss rate and low wind velocity, the Fe-K line centroid 
approaches the collisional ionization equilibrium value in 
several hundred years While the measured line centroids 
as a function of SNR radius do generally agree with the 
models, the observed line luminosities for as many as half 
of the sampled remnants are greater than the modeled 
luminosities by as much as an order of magnitude. We 
discuss these results in more detail below. 


et al. (2013) postulated that W49B is the result of a bipo¬ 
lar Type Ib/c SN, based on the morphology and other 
spectral characteristics. However, the low CSM densities 
inferred from X-ray and radio obs ervat i ons to be a round 
Ib/c progenitors (e.g., SN 2007gr; ISoderberg et al.|20To| 
appears incongruent with a high Fe-K charge state and 
line luminosity in the spectrum of W49B. Type Ib/c SNe 
probably result from Wolf-Rayet stars with fast winds 
and high mechanical luminosities. Compared to red and 
yellow supergiant progenitors, the 100x faster wind ve¬ 
locities in WR progenitors transport much of the cir- 
cumstellar material to larger distances, leaving behind 
a largely low density CSM cavity, though evidence does 
exist for more complex environments around some SN Ic, 
such as in the case of SN 2007bg, which appears to have 
und erwent sever al differing mass- loss phases prior to the 


SN (Salas et al. 


zawa et al. 


(2009) found evidence for an 


Additionally, 

overionized plasma in W49 B, as evidenced by the ratio of 
H- to He-li ke lines (see also Lopez et al. 2013 Yamaguchi] 


Chevalier (20051 typed several Galactic SNRs based 
on the properties of the swept up circumstellar mate¬ 
rial, classifying many as either Type Ilb/L or Type IIP. 
Main sequence progenitors for Type IIP and Ilb/L SN 
span a mass range of 10-25 M 0 , and have slow dense 
winds resulting in much of the lost mass remaining close 


et al. 
fdTS] 


l. 2014]). 
TNRs th. 


Moriya (20121 recently presented a model 


that show overionized plasmas, and suggested 


that red supergiants can deposit the mass required for 
rapid iouizatioir close enough to the progenitor to ove¬ 
rionize the shocked material at early times. He noted 
that the time required to reach ionization equilibrium in 
a RSG wind is ~ 100x faster than in a W-R wind, un- 
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less collisional ionization equilibrium in the W-R wind 


occur s almost immediately after the explosion (Moriya 


2012). However, an emergent class of supernovae typed 


as I bn has recently been discovered that displays a strong 
circumstellar interaction like those found in Iln’s as well 
as no H in th e early time o ptical spectra, li ke lb’s, such 
as SN 2006cj ( Chugai|[2009 ) and SN 2011w (Smith et al. 


2012|. The progenitors of these systems include Wolf- 


Rayet and LB Vs with zero age main sequence masses 
M > 40 M 0 , a nd pos sibly display eruptive events just 
prior to the SN 
may explain bo 


Smith (2014). A progenitor such as this 
h the overionization and high Fe abun¬ 
dance observed in W49B and reconcile these apparently 
contradictory results with theory. In any event, detailed 
modeling of a W49B-like progenitor, explosion, and sub¬ 
sequent evolution might be required to explain this odd 
object. 


4.2. Spectral vs. Dynamical Quantities 

As seen in Figure^ (right), our modeled line luminosi¬ 
ties fall below the observed values for many remnants. 
Here we argue that this may be associated with issues 
related to the very end stages of the progenitor’s evolu¬ 
tion. While the observed SNR radii appear consistent 
with the modeled radii, the luminosities differ by more 
than an order of magnitude. So the question is whether 
one can increase the luminosity without strongly impact¬ 
ing the radius. For core-collaps e super nov ae, t he forward 
shock radius R b cx M _1 /(«- s ) ( |Chevalier|| 19821 ). Thus if 
one were to increase the mass-loss rate prior to the super¬ 
nova, we can estimate the change in bl astwave radius as 
a function of change in M. First, from 
the blastwave radius is: 


Chevalier 


(1982), 


R b oc 


Ag n 

Q 


l/(n-s) 


71-3 

t”- S , 


( 2 ) 


where A is a constant, dependent upon the shape of the 
ejecta and circumstellar environment, n is the power law 
index that describes the shape of the ejecta, and s is the 
power law index that describes the circumstellar environ¬ 
ment. In this paper, n = 9 and s = 2. g n is a constant, 
dependent upon the ejecta mass and explosion energy, q 
is defined as M /4ttv w . If the mass-loss rate is increased 
by some amount x prior to the supernova ( q' = \g) , such 
that the blastwave radius is now R' b oc gl -1 /!"- 8 ))', then 
the fractional change in the radius, over the extent of 
time that the blastwave moves through the circumstellar 
shell is: 


R b -R' b _ g _1 A n - s ) _ g(—!/("—*))' 

R b ~ q-D(n-s) W 

= 1 - x~ 1/(n ~ s) ■ (4) 

For large changes in the mass-loss rate, this effect can 
be quite large. For instance, an order of magnitude in¬ 
crease in the mass-loss rate results in a 30% decrease in 
the blastwave radius. However, we do not expect this in¬ 
creased mass-loss to persist for very long, so the extent 
of the circumstellar medium occupied by the increased 
density is small, essentially amounting to a thin circum¬ 
stellar shell. After the blastwave breaks through this 


shell, it will accelerate into the lower density slow RSG 
wind. 

In contrast, the X-ray luminosity Lx ~ n^A (T)V, 
where n e is the number of free electrons, A(T) is the 
cooling function, and V is the emitting volume. Since 
n e oc M, Lx ~ M 1 . A three-fold increase in the mass- 
loss rate will result in an approximately order of magni¬ 
tude increase in L y , an increase in collisional ionization, 
and thus an increase in Fe-K line emission. Thus, the 
problem posed by Figure [2] (right) can be qualitatively 
addressed by altering the density of the circumstellar en¬ 
vironment close to the progenitor. This could be viewed 
as akin to higher mass-loss just prior to the supernova. 
We note that this argument may break down for shells 
that are sufficiently dense to cause significant radiative 


losses (Crowther 2007 Pan et al.||2013 Smith||2014|). 

In Figure[2j(left) we plot the modeled and observed Fe- 
K line centroid versus SNR radii. As shown in Figure [TJ 
the energy centroids for both the model and data are 
generally greater than 6550 eV. In general, the observed 
and modeled line centroids as a function of energy are 
in agreement with one another. However, some SNRs, 
G292.0+1.8 in particular, do not agree with any of the 
models. Interestingly, in the line luminosity versus radii 
plots, G292.0+1.8 also does not agree with the models. 
This may point to a more complicated circumstellar en¬ 
vironment around G29 2.0+1.8, evidence of which is seen 
in X-ray observations ( Park et al.|2002 Lee et al.|2010 ). 
However, in terms of the models, we find that they gen- 
erally predict Fe-K line centroids in excess of 6550 eV 
for a broad range of radii. For the high mass-loss rate 
models, the Fe-K line centroid is approaching collisional 
ionization equilibrium - increasing the mass-loss rate to 
affect a change in the luminosity would not alter this 
conclusion. 


5. CONCLUSIONS 

We have presented an initial attempt to connect some 
properties of core-collapse supernova, namely the com¬ 
position and a parameterization of the circumstellar en¬ 
vironment, to some observable bulk properties of super¬ 
nova remnants. We find that: 

1. Using a model for the dynamical and spectral evo¬ 
lution of shocked supernova ejecta and swept up 
material, we computed Fe-K line centroids and lu¬ 
minosities for a broad range of parameters, includ¬ 
ing explosion energetics, ejecta mass, and circum¬ 
stellar environment. We found that they distinctly 
differ from those same properties computed from 
models for Type la supernovae. 

2. We compared our models to measured Fe-K line 
properties and found general agreement between 
the models and observations. For those SNRs that 
do not show good agreement with these simple 
models we put forth possible reasons for the dis¬ 
crepancy. 

3. We discuss the relationship between SNR radius 
and Fe-K luminosity, and propose that higher 
mass-loss rates prior to the supernova may be 
able to reconcile the discrepancy between the mod¬ 
eled and observed luminosities. The required 
































6 


Patnaude et al. 


higher mass-loss may arise either through erup¬ 
tive, episodic events prior to the SNe, such as in 
SN 2009ip, or increased clumping in the wind dur¬ 
ing the later stages of evolution. Both of these 
possibilities would result in a stellar wind profile 
that differs from the r~ 2 winds considered here. 

For core-collapse supernovae, the range of ejecta mass, 
explosion energetics, and circumstellar environment can 
be large. However, canonical values for supernova pa¬ 
rameters can reproduce many observations. Objects such 
as Cas A and W49B require more detailed models, and 
any other SNR shown in Figure [l] may be scrutinized in 
more detail in order to pin down the exact parameter 
space. In terms of connecting supernovae (and their pro¬ 
genitor) models to supernova remnants, more detailed 
studies of young SNRs such as SN 1993J, SN 1996c, and 
NGC 4449-1 are required in order to bridge the gap be¬ 
tween truly young objects and Galactic remnants. Ad¬ 
ditionally, studies of these objects will allow us to probe 
the mass-loss history of the progenitor at times < 10 5 
years before the explosion. Future concept X-ray mis¬ 
sions with both high imaging and spectral resolution, 
such as Smart- fl will foster the detection and study of 
nearby extragalactic SNRs which can be added to the 
sample discussed here. 
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Fig. 1.— Left: Fe-K line luminosity versus centroid energy for Galactic and Magallenic Cloud la and core-collapse supernova remnants. 
The transparent shaded regions correspond to the models listed in Table[l]with yellow corresponding to model D DTa, light b lue to DDTg, 
red to 1987A, green to 1993J, magenta to sl2D and dark blue to s25D. The data are taken from Table 1 of|Yamaguchi et al.112014|. Right: 
Fe-K luminosity as a function of scaled SNR age. The solid curves are the FS emission, while the dashed curves are tor tne RS emission. 
The two sets of curves represent the two sets of models listed in Table [T] The SN 1987A model does not produce any appreciable Fe-K 
emission from the RS. 


TABLE 1 

Model Parameters 


Ejecta Model 

Esiv 

10 s1 erg 

M ej 
M 0 

n amb 

_ Q 

cm 

Vwind 
km S —1 

M b 

10 -5 Mq yr -1 

Ref. 


DDTa 
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1.38 

0.1 - 3.0 



Badenes et al. 

(2008 


DDTg 
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1.38 
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10-20 

1-2 

1-2 

Badenes et al. 

(2008 


sl2D 

s25D 

1.21 

1.21 

8.87 

12.2 


ih 

Thi 

s worK 
s work 

1987A 

1.10 

14.7 


10-20 

1-2 

|Saio et al. 

(1988 


1993J 

2.00 

2.92 


10-20 

1-2 

Nozawa et al. 

(2010 



a Type la models evolve in a constant density environment. 

b Core—collapse models evolve in a circumstellar environment shaped by isotropic 
mass-loss. 
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Fig. 2. — Left : Fe-K line centroid energy as a function of radius for core-collapse supernova remnant measurements and CCSNe models 
listed in Table |T[ Right : Fe-K line luminosity as a function of radius for core-collapse supernova remnant measurements and CCSNe models 
listed in Tabled 







































